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<: 

Microscopic black holes explode with their temperature varying inversely 
fN| . as their mass. Such explosions would lead to the highest temperatures in the 

present universe, all the way to the Planck energy. The possibility that a 
quasi-stationary shell of hot matter surrounds these black holes has recently 

■ been proposed and studied with relativistic Boltzmann transport equations 
and with relativistic viscous fluid dynamics. For example, a black hole with a 

' mass of 10 10 g has a Hawking temperature of 1 TeV, a Schwarszchild radius 

■ of 1.6xl0 -5 fm, a luminosity of 7xl0 27 erg/s, and has less than 8 minutes 
^ c*j to live. It is an outstanding theoretical challenge to describe the conditions 

exterior to such microscopic black holes and a great challenge to finally detect 
Q ■ them in the new millennium. 

a: 

Hawking radiation from black holes Jl| is of fundamental interest because it 
relies on the application of relativistic quantum field theory in the presence of the 
strong field limit of gravity, a so far unique situation. It is also of great interest 
because of the temperatures involved. A black hole with mass M radiates thermally 
with a Hawking temperature 



mp 



where mp = G^ 1 ^ 2 = 1.22 x 10 19 GeV is the Planck mass. (Units are h = c = 
/cb = 1.) In order for the black hole to evaporate it must have a temperature 
greater than that of the present-day black-body radiation of the universe of 2.7 
K = 2.3 xlO~ 4 eV. This implies that M must be less than 1% of the mass of the 
Earth, hence the black hole most likely would have been formed primordially and 
not from stellar collapse. The black hole temperature eventually goes to infinity 
as its mass goes to zero, although once Th becomes comparable to the Planck 
mass the semi-classical calculation breaks down and the regime of full quantum 
gravity is entered. Only in two other situations are such enormous temperatures 
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achievable: in the early universe (T similarly asymptotically high) and in central 
collisions of heavy nuclei like gold or lead (T = 500 MeV is expected at the RHIC 
(Relativistic Heavy Ion Collider) just completed at Brookhaven National Laboratory 
and T = 1 GeV is expected at the LHC (Large Hadron Collider) at CERN to be 
completed in 2005). Supernovae and newly formed neutron stars are unlikely to 
ever exceed a temperature of 50 MeV. To set the scale from fundamental physics, 
we note that the spontaneously broken chiral symmetry of QCD gets restored in 
a phase transition/rapid crossover at a temperature around 160 MeV, while the 
spontaneously broken gauge symmetry in the electroweak sector of the standard 
model gets restored in a phase transition/rapid crossover at a temperature around 
100 GeV. The fact that temperatures of the latter order of magnitude will never 
be achieved in a terrestrial experiment should motivate us to study the fate of 
primordial black holes during the final minutes of their lives when their temperatures 
have risen to 100 GeV and above. The fact that primordial black holes have not yet 
been observed should not be viewed as a deterrent, but rather as a challenge! 

There are at least two intuitive ways to think about Hawking radiation from 
black holes. One way is vacuum polarization. Particle-antiparticle pairs are con- 
tinually popping in and out of the vacuum, usually with no observable effect. In 
the presence of matter, however, their effects can be observed. This is the origin of 
the Lamb effect first measured in atomic hydrogen in 1947. When pairs pop out of 
the vacuum near the event horizon of a black hole one of them may be captured by 
the black hole and the other by necessity of conservation laws will escape to infinity 
with positive energy. The black hole therefore has lost energy - it radiates. Due to 
the general principles of thermodynamics applied to black holes it is quite natural 
that it should radiate thermally. An intuitive argument that is more quantitative 
is based on the uncertainty principle. Suppose that we wish to confine a massless 
particle to the vicinity of a black hole. Given that the average momentum of a 
massless particle at temperature T is approximately 7rT, the uncertainty principle 
requires that confinement to a region the size of the Schwarzschild diameter places 
a restriction on the minimum value of the temperature. 

7rT-2r 5 >l/2 (2) 

The minimum is actually attained for the Hawking temperature. The various phys- 
ical quantities are related as = 2M/rrip = l/AirT^. 

The number of particles of spin s emitted with energy E per unit time is given 
by the formula 

dN s _ r, 1 

dEdt 2tt exp(E/T H ) - (-l) 2s ' [ ' 

All the computational effort really goes into calculating the absorption coefficient 
T s from a relativistic wave equation in the presence of a black hole. Integrating over 
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all particle species yields the luminosity. 



L 



dM 



a(M) 




U-K 2 a{T H )T 2 H 



(4) 



dt 



M 2 



Here a(M) is a function reflecting the species of particles available for creation in 
the gravitational field of the black hole. It is generally sufficient to consider only 
those particles with mass less than T#; more massive particles are exponentially 
suppressed by the Boltzmann factor. Then 



where N s is the net number of polarization degrees of freedom for all particles 
with spin s and with mass less than Th- The coefficients for spin 1/2, 1 and 2 
were computed by Page || and for spin by Sanchez J4|. In the standard model 
N = 4 (Higgs), Ni/2 = 90 (three generations of quarks and leptons), N% = 24 
(SU(3) xSU(2) xU(l) gauge theory), and N 2 = 2 (gravitons). This assumes T H is 
greater than the temperature for the electroweak gauge symmetry restoration. Nu- 
merically a(Tn > 100 GeV) = 4.43 x 10~ 3 . Starting with a black hole of temperature 
T H , the time it takes to evaporate/explode is 



This is also the characteristic time scale for the rate of change of the luminosity of 
a black hole with temperature Th- 

At present a black hole will explode if Th > 2.7 K. This corresponds to a critical 
mass of M < 4.6 x 10 25 g which is approximately 1% of the mass of the Earth. More 
massive black holes are cooler and therefore will absorb more matter and radiation 
than they radiate, hence grow with time. Taking into account emission of gravitons, 
photons, and neutrinos a critical mass black hole today has a Schwarszchild radius 
of 68 microns and a lifetime of 2 x 10 43 years. 

There is some uncertainty over whether the particles scatter from each other 
after being emitted, perhaps even enough to allow a fluid description of the wind 
coming from the black hole. Let us examine what might happen as the black hole 
mass decreases and the associated Hawking temperature increases. 

When Th <C m e (electron mass) only photons, gravitons, and neutrinos will 
be created with any significant probability. These particles will not interact with 
each other but will be emitted into the surrounding space with the speed of light 
[p], [|. Even when Th ~ Tn e the Thomson cross section is too small to allow the 
photons to scatter very frequently in the rarified electron-positron plasma around 
the black hole. This may change when Th ~ 100 MeV when muons and charged 
pions are created in abundance. At somewhat higher temperatures hadrons are 
copiously produced and local thermal equilibrium may be achieved, although exactly 



a = 2.011 x 10~ s 4200iVo + 2035iVi /2 + 835iVi + 95iV 2 



(5) 



At 




(6) 



3a(T H )(87rT H ) 3 ' 
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how is an unsettled issue. Are hadrons emitted directly by the black hole? If so, 
they will be quite abundant at temperatures of order 150 MeV because their mass 
spectrum rises exponentially (Hagedorn growth as seen in the Particle Data Tables 
13]). Because they are so massive they move nonrelativistically and may form a very 
dense equilibrated gas around the black hole. But hadrons are composites of quarks 
and gluons, so perhaps quarks and gluon jets are emitted instead? These jets must 
decay into the observable hadrons on a typical length scale of 1 fm and a typical 
time scale of 1 fm/c. This was first studied by MacGibbon and Webber JJJ and 
MacGibbon and Carr ][7). Subsequently Heckler || argued that since the emitted 
quarks and gluons are so densely packed outside the event horizon that they are not 
actually fragmenting into hadrons in vacuum but in something more like a quark- 
gluon plasma, so perhaps they thermalize. He also argued that QED bremsstrahlung 
and pair production were sufficient to lead to a thermalized QED plasma when Th 
exceeded 45 GeV ||. These results are somewhat controversial and need to be 
confirmed. The issue really is how to describe the emission of wavepackets via the 
Hawking mechanism when the emitted particles are (potentially) close enough to 
be mutually interacting. A more quantitative treatment of the particle interactions 



on a semiclassical level was carried out by Cline, Mostoslavsky and Servant flO 
They solved the relativistic Boltzmann equation with QCD and QED interactions in 
the relaxation-time approximation. It was found that significant particle scattering 
would lead to a photosphere though not perfect fluid flow. 

The picture proposed by Heckler is as follows: Hawking radiation is emitted as 
usual in the vicinity of the Schwarzschild radius. Bremsstrahlung scattering leads 
to approximate thermalization and fluid flow at a radius r > r$- Thermal contact 
is lost at a photosphere radius of r p > r . By then most electrons and positrons 
have annihilated; the energy goes off to infinity in the form of photons with average 
energy much less than Th- 

Figures 1-3 show some results from Cline et al. The number of scatterings in 
the QCD shell of matter interior to the hadronization radius grows significantly 
with increasing Th- Indeed the total increase in particle production due to QCD 
bremsstrahlung and to quark and gluon fragmentation goes linearly with Th, while 
the outer radius of the QCD shell of matter grows logarithmically. 

Rather than pursuing the Boltzmann transport equation I will now attempt a 
description in terms of relativistic viscous fluid dynamics |TT| . Let us assume that a 



primordial black hole is surrounded by a shell of expanding matter in approximately 
local thermal equilibrium when Th is large enough. A detailed description of how 
this situation comes to be is a difficult problem as discussed above and needs further 
investigation. The applicability can and will be checked a posteriori. The relativistic 
imperfect fluid equations describing a steady-state, spherically symmetric flow with 
no net baryon number or electric charge and neglecting gravity (see below) are 
T^ v . v = black hole source. The nonvanishing components of the energy-momentum 
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Figure 1: Average number of scatterings in QCD photosphere as a function of radius 
for several black hole temperatures. From [PTj| . 
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Figure 2: Radius of the outer surface of the QCD photosphere versus logarithm of 



the black hole temperature. From |L0 . 
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Figure 3: Total particle production in the QCD photosphere versus black hole tem- 
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tensor in radial coordinates are 12 



rpOO 
rpOr 

rprr 



1 2 (P + e)-P + v 2 AT diss 
v 7 2 (P + e)+vAT diss 
v 2 ^(P + e) + P + AT ( 



diss 



(7) 



representing energy density, radial energy flux, and radial momentum flux, respec- 
tively, in the rest frame of the black hole. Here v is the radial velocity with 7 the 
corresponding Lorentz factor, u = wy, e and P are the local energy density and 
pressure, and 

4 / du u\ „ 9 / du 2u\ 



AT, 



diss 



-^7 



dr 



dr 



+ 



(8) 



where 77 is the shear viscosity and ( is the bulk viscosity. A thermodynamic identity 
gives Ts = P + e for zero chemical potentials, where T is temperature and s is 
entropy density. There are two independent differential equations of motion to solve 
for the functions T(r) and v(r). These may succinctly be written as 



Tr ^) 








(9) 
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An integral form of these equations is probably more useful since it can readily 
incorporate the input luminosity Li from the black hole. The first represents the 
equality of the energy flux passing through a sphere of radius r with the luminosity 
of the black hole. 

47rr 2 T 0r = ^ ( 1Q ) 

The second follows from integrating a linear combination of the differential equa- 
tions. It represents the combined effects of the entropy from the black hole together 
with the increase of entropy due to viscosity. 



4nr 2 us = 4ir I dr' r' 2 



T 




+ ^ (11) 



T> 



The term Li/Tn arises from equating the entropy per unit time lost by the black hole 
—dSbh/dt with that flowing into the matter. Using the area formula for the entropy 
of a black hole, Sbh = mp7rr| = 4irM 2 /rrip, and identifying —dM/dt with the lumi- 
nosity, the entropy input from the black hole is obtained. The above pair of equations 
are to be applied beginning at some radius r« greater than the Schwarzschild radius 
rs, that is, outside the quantum particle production region of the black hole. The 
radius r\ at which the imperfect fluid equations are first applied should be chosen 
to be greater than the Schwarzschild radius, otherwise the computation of particle 
creation by the black hole would be invalid. It should not be too much greater, 
otherwise particle collisions would create more entropy than is accounted for by the 
equation above. The energy and entropy flux into the fluid come from quantum 
particle creation by the black hole at temperature T#. Gravitational effects are of 
order rs/r, hence negligible for r > (5 — 10) rs. 

Determination of the equation of state as well as the two viscosities for tem- 
peratures ranging from MeV to TeV and more is a formidable task. Here we shall 
consider two interesting limits and then a semi-realistic situation. A realistic, quan- 
titative description of the relativistic black hole wind, including the asymptotic 
observed particle spectra, is currently under investigation by my graduate student 
Ramin Daghigh and myself. 

First, consider the nonviscous limit (like milk) with an equation of state e = aT 4 , 
s = (4/3) aT 3 , and rj = ( = 0. This is equivalent to assuming that the mean free 
paths of the particles are all small compared to the length scale over which the 
temperature and other thermodynamic quantities change significantly. A scaling 
solution, valid when 7 3> 1, is T(r) = T (r /r) and 7(7") = 7o(r/r ), where 7oT = 
Tu- The tq is any reference radius satisfying the stated criterion. 

Second, consider the highly viscous isoergic limit (like honey) in the sense that 
the flow velocity approaches a limiting value v at large r. This requires a power- 
like equation of state e cx T 5 and viscosities rj oc ( oc T & / 2 . It results in the 
scaling solution T(r) = T (r /r) 2 / 5 . This is not very realistic: a massless gas with 
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dimensionless coupling constants and 5 = 4 would require viscosities of order T 2 
whereas one would expect T 3 on dimensional grounds. 

Now consider a semi-realistic situation with e = aT 4 , s = (4/3)aT 3 , rj = bsT 3 , 
and ( = bsT 3 . This is typical of relativistic gases with dimensionless coupling 
constants, although quantum effects will give logarithmic corrections [L4|. A 
scaling solution, valid at large radii when 7 ^> 1, is T(r) = T (r /r) 2 / 3 and 7(7") = 
7o(r/r ) 1 / 3 . This r-dependence of T and 7 is exactly what was conjectured by 
Heckler §. 

Is the semi-realistic situation described above really possible? Can approximate 
local thermal equilibrium, if once achieved, be maintained? The requirement is that 
the inverse of the local volume expansion rate 9 = u^.^ be comparable to or greater 
than the relaxation time for thermal equilibrium jl2| . Expressed in terms of a local 
volume element V and proper time r it is 9 = (l/V)dV/d,T, whereas in the rest frame 
of the black hole the same quantity can be expressed as (l/r 2 )d(r 2 u)/dr. Explicitly 

e = 7 -IiC<±) 2l3 = ^LT. (12) 

3r V r / 3r T 

Of prime importance in achieving and maintaining local thermal equilibrium in a 
relativistic plasma are multi-body processes such as 2 —> 3 and 3 — ► 2, etc. This has 
been well-known when calculating quark-gluon plasma formation and evolution in 
high energy heavy ion collisions |15[ and has been emphasized in ref. |8], [| in the 
context of black hole evaporation. This is a formidable task in the standard model 
with its 16 species of particles. Instead we make three estimates for the requirement 
that local thermal equilibrium be maintained. The first and simplest estimate is to 
require that the thermal DeBroglie wavelength of a massless particle, 1/3T, be less 
than 1/9. The second estimate is to require that the Debye screening length for each 
of the gauge groups in the standard model be less than 1/9. The Debye screening 
length is the inverse of the Debye screening mass where n = 1,2,3 for the gauge 
groups U(l), SU(2), SU(3). Generically m° oc g n T where g n is the gauge coupling 
constant and the coefficient of proportionality is essentially the square root of the 



number of charge carriers fIB| . For example, for color SU(3) = gsyjl + N[/6T 
where iVf is the number of light quark flavors at the temperature T. The numerical 
values of the gauge couplings are: g\ = 0.344, g 2 = 0.637, and #3 = 1.18 (evaluated 
at the scale rriz) f|. So within a factor of about 2 we have m D « T. The third 
and most relevant estimate is the mean time between two-body collisions in the 
standard model for temperatures greater than the electroweak symmetry restoration 
temperature. This mean time was calculated by Carrington and me [0 in the 
process of calculating the viscosity in the relaxation time approximation. Averaged 
over all particle species in the standard model one may infer from that paper an 
average time of 3.7/T. Taking into account multi-body reactions would decrease that 
by about a factor of two to four. All three of these estimates are consistent within 
a factor of 2 or 3. The conclusion to be drawn is that local thermal equilibrium 
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Figure 4: Entropy density as a function of temperature, excluding neutrinos and 
gravitons. It is assumed that the QCD phase transition is first order and the EW 
phase transition is second order. 

should be achieved when 9 < T. Once thermal equilibrium is achieved it is not lost 
because 9/T is independent of r. The picture that emerges is that of an imperfect 
fluid just marginally kept in local equilibrium by viscous forces. 

The hot shell of matter surrounding a primordial black hole provides a theoretical 
testing ground rivaled only by the big bang itself. In addition to the questions 
already raised, one may contemplate baryon number violation at high temperature 
and how physics beyond the standard model might be important in the last few 
minutes in the life of a primordial black hole. To illustrate this I have plotted a semi- 
realistic parametrization of the equation of state in figure 4. Gravitons and neutrinos 
are not included. I have assumed a second order electroweak phase transition at a 
temperature of 100 GeV. Above that temperature the standard model has 101.5 
effective massless bosonic degrees of freedom (as usual fermions count as 7/8 of a 
boson). There is a first order QCD phase transition at a temperature of 160 MeV. 
The number of effective massless bosonic degrees of freedom changes from 47.5 just 
above this critical temperature (u, d, s quarks and gluons) to 7.5 just below it 
(representing the effects of all the massive hadrons in the particle data tables) [17] . 
Below 30 MeV only electrons, positrons, and photons remain, and finally below a 
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few hundred keV only photons survive in any appreciable number. This is a rather 
crude description of the real equation of state that is still a matter of theoretical 
uncertainty. One needs knowledge of s(T),r)(T)X(T) over a huge range of T! Of 
course, these are some of the quantities one hopes to obtain experimental information 
on from observations of exploding black holes. 

The parallelism with high energy nucleus-nucleus collisions cannot go unnoticed 
and without a certain amount of cross fertilization. In heavy ion physics the best 
theories are founded on relativistic fluid dynamics and on microscopic transport 



models fll8| . Both approaches can be fruitfully employed to describe the conditions 
outside the event horizon of a microscopic black hole. There is much work to be done 
in this area. Finally, such black holes may contribute to the highest energy cosmic 
rays whose origin is a long-standing puzzle. Experimental discovery of exploding 
black holes will be one of the great challenges in the new millennium! 
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